Branching ratios of radiative transitions in O VI 
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Abstract 

We study the branching ratios of the allowed and forbidden radiative transitions among the first few (9) fine structure 
levels of O VI using relativistic coupled cluster theory. We find irregular patterns for a number of transitions with in 
n-complexes with n < 4. We have used the exisiting values of the allowed electric dipole {El) transition as a benchmark 
of our theory. Good agreement with the existing values establish accuracies of not only the theoretical method but the 
basis function as well. In general the electric quadrupole {E2) transition probabilities are greater in magnitude than 
magnetic dipole (Ml) transition probabilities, whereas for medium atomic transition frequencies they are of the same 
order of magnitude. On the other hand if the transitions involved are between two fine structure components of the same 
term, then the Ml transition probability is more probable than that of E2. We have analyzed these trends with physical 
arguments and order of magnitude estimations. The results presented here in tabular and graphical forms are compared 
with the available theoretical and observed data. Graphical analysis helps to understand the trends of electric and magnetic 
transitions for the decay channels presented here. Our calculated values of the lifetimes of the excited states are in very 
good agreement with the available results. 

PACS numbers : 31.15.Dv, 31.10.+Z, 32.70.Cs 
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\0 % Introduction 
00 . 

Comprehensive and accurate transition probabiHty data sets are needed to determine the abundances of neutral oxygen 
and all of its ions in different astrophysical spectra. Determination of radiative lifetimes can provide the absolute scale 
C ■ for converting the branching fractions into atomic transition probabilities and vice-versa. The emission line diagnostics of 
allowed and forbidden transitions for VI have been used for a wide range of spectral modeling of planetary nebulae and H 
11 regions [U and in NGC 2867 [3]. The range of these spectral studies lies between ultraviolet (UV) and optical region. 
. . The most prominent emission Hue doublets of O VI are (AA 1032, 1038) and (AA 3811, 3834). They lie in the far-UV and 
. ^ Optical regions of the spectra respectively in many central stars. Moreover O VI is one of the highest observed states of 
^ ionization in stellar spectra and is thus considered a rich source of data for high temperature plasma diagnostics. Tayal [4j 
$—1 and Aggarwal et al [5j provide more references. Tayal [4] has studied the allowed transitions and electron impact collision 
. 5^ strengths in the Breit-Pauli R-matrix framework whereas Aggarwal et al. [5] have employed the fully relativistic GRASP 
[6] and DARC ^ codes for the same study. Froese Fischer et al. have employed multi-configuration Hartree-Fock (MCHF) 
method to calculate the allowed (El) transition probabilities and the lifetime of the excited states [8]. 

The calculation of radiative transitions for many-electron atoms requires a good description of the electron correlation 
which determines the movement of the radiating electrons. In this article we have used relativistic open-shell coupled cluster 
method (RCCM), one of the most accurate theories to describe the electron correlation effects in many-electron atoms, to 
calculate allowed and forbidden radiative transition probabilities. 

This article deals with the study of allowed and forbidden transitions in O VI ion using relativistic coupled cluster (RCC) 
theory. The present study of the allowed and forbidden transitions is used to identify the dependencies of the line strengths, 
transition probabilities and oscillator strengths on the transition energies. From the definition, the dependencies of transition 
probabilities and oscillator strengths on the transition energies is quite clear but we have observed that the line strengths also 
follow a trend which can be understood by studying the physics involved in the problem. Graphical analysis of the results 
presented in this article can be used to understand the strengths of different electromagnetic transitions for a particular 
emission. 
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2 Theoretical methods 



In this work we have employed RCCM to study the branching ratios of VI. The virtue of RCCM is that it is an all- 
order non-perturbative theory and size extensive in nature [9]. Higher order electron correlation effects can be incorporated 
more efficiently than by using order-by-order diagrammatic many-body perturbation theory (MBPT) . Extensive literature is 
available on RCC theory and application to calculate atomic properties [IHlIll], so it will not be provided here. Our method 
of calculation is described in details in a recent review article [12] . 



2.1 Radiative transitions 

The interaction of electromagnetic radiation with atoms in most plasma applications is generally treated in what is known 
as first quantization, i.e. the quantized energy levels of the atom interact with a continuous radiation field. For spontaneous 
decay, transition probabilities (A) of different electromagnetic multipoles are well known. Oscillator strengths (f ) and hence 
the weighted oscillator strengths (gf ), which are the products of the degeneracies of the initial states of the transition and f 
are very important in astrophysics and can be used to obtain information about the spectral formation and other quantities 
of astrophysical interests. The knowledge of transition probabilities is often used to determine the line ratios of the density 
diagnostics of the plasmas. 

In this article Ek{Mk) stands for electric (magnetic) type transitions and k — 1,2.. correspond to dipole, quadrupole 
transition respectively. Here s^'^Af^ are the line strengths; A is the wavelength in A and t^^/ is the statistical weight for the 
upper level for a transition / — > i. Afi and f/i are the transition probabilities and oscillator strengths respectively. The A 
and f values for different types of transitions with the line strength (S) are related by the following standard equations : 
for the electric dipole (El) transitions : 

2.0261 X IQis „ 303.75 ei 

for the electric quadrupole {E2) transitions : 

1.1199 X IQis 167.89 
A/^ = -T5 ^ and f = ^^S^^ (2) 



UJ. 



for the magnetic dipole (Ml) transitions : 

A/» = ^'""',',3" S^^'^ and f/, = ----- - gMi (3-) 



2.6974 x 10^3^^^^^^^ 4.044,10-3^^^^ 



and for the magnetic quadrupole (M2) transitions : 

_ 1.4910x10^3 2.236 x 10-3 ^,,, 

- cofX% ^ - X%u, ^ ■ 

The single particle reduced matrix elements for the El, E2, Ml and M2 transitions are given by, 

{Kf\\el\\K,) = IIj.) x I r{PfP, + QfQi)dr (5) 

{KfW e2 11^,) = i| {jfW ) IIj-) X / j2(fcr) (P/P, + QjQ,) dr, (6) 



(K/||ml||«:.) = A(j^l|CWl|j-)x j 11±^^,[kr){PfQ., + QfP,)dr (7) 

and 

{nf\\m2\\n,)^^{jf\\Cf^\\j,)x j ^l±^^,ikr) {PfQ, + QjP,) dr (8) 

respectively. Here j's and k's are the total orbital angular momentum and the relativistic angular momentum quantum 
numbers respectively; k is defined as Loa where ui is the single particle difference energy and a is the fine structure constant. 
The single particle orbitals are expressed in terms of the Dirac spinors with Pi and Qi as the large and small components for 
the ith spinor, respectively. The angular coefficients are the reduced matrix elements of the spherical tensor of rank m and 
are expressed as 
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Table 1: Excitation Energies for different states (in Ryd) for VI. RCCM stands for the present calculation. GRASP and 
CIVS are the calculations done by Aggarwal et al [5j and Tayal [4] respectively. The numbers presented under the column 
NIST are obtained from the database |18l. 



Index States 



RCCM NIST GRASP CIV3 



ls23p2p3/2 
1523^21)3/2 

8 ls23d2D5/2 
9 



0.0000 
0.8790 
0.8845 
5.8337 
6.0735 
6.0751 
6.1502 
6.1506 
7.7732 



0.0000 
0.8782 
0.8831 
5.8325 
6.0701 
6.0715 
6.1476 
6.1481 
7.7703 



0.0000 
0.8863 
0.8910 
5.8251 
6.0646 
6.0660 
6.1391 
6.1395 
7.7615 



0.0000 
0.8876 
0.8944 
5.8292 
6.0676 
6.0696 
6.1421 
6.1427 
7.7663 



i^fW ll'^.) = (-1^+^/^(2.7/ + l)(2j. + 1) ( I ^\y{lf,m,k), (9) 

with 

J 1 if + m + even 

otherwise ^1°) 

and I being the orbital angular momentum quantum number. When z — kr is sufHciently small the spherical Bessel function 
j„(z) is approximated as 

^"^^^"^ (2n+l)!! " l-3-5----(2n+l)- ^^^^ 



3 Results and discussions 

In this manuscript we have studied different electromagnetic transitions in VI. We have applied RCCM to calculate the 
energy levels, Hne strengths, transition probabilities and oscillator strengths for the first nine levels. The calculation of 
radiative transitions in many-electron atoms requires a good numerical description of the electronic states and the wave 
functions. From the excitation energies presented in table [Hit is evident that determinations of the energy levels are very 
accurate with an average error of less than 0.2% which is presented graphically in figure[2l The comparison with the presently 
available theoretical and observed data establishes the accuracy of RCCM to determine the state energies which are used to 
calculate the transition properties related to astrophysical interests. 

In the actual computation, the Dirac-Fock (DF) ground state and excited state properties of O VI are computed using the 
finite basis set expansion method (FBSE) ^3j of a large basis set of (35s32p25d25/) Gaussian functions with the assumption 
that the nucleus has a Fermi type finite structure. Excitations from all core orbitals have been considered within the 
active space consisting of 14sl3pll(il0/ orbitals. Table [2] contains the comparative data for the allowed electric dipole (El) 
transitions for all the nine levels. We have emphasized the two most important fine structure (FS) doublets ((AA1032, 1038) 
and (AA 3811, 3834)). The (AA1032, 1038) FS doublet is the strongest emission feature in the far-UV spectra of RR Tel [H] 
and from the Berkeley spectrograph on ORFEUS [H]. They are also observed in many other symbiotic stars [4]. On the 
other hand the (AA 3811, 3834) FS doublet Hes in the optical region of many central stars and the presence of this feature 
is often referred to as "0 VI sequence". A close comparison of the allowed transition (El) properties with the GRASP and 
CIV3 calculations shows that this is the most accurate determination of the allowed transition probabilities in VI to date. 

Tables [HM] contain the details of the forbidden transitions in O VI. This is the first theoretical determination of forbidden 
transition probabilities and oscillator strengths for O VI. They are related to the corresponding Hne strengths by Eqs. (I2][4]). 
From the data presented in tables |[2][5]) we have determined the radiative lifetimes of the excited states which decay via 
allowed and forbidden transitions to lower levels. We have observed that among the forbidden transitions between the three 
fine structure doublets, the E2 line-strengths {Se2) are greater than Ml line-strengths (Smi) except for ls^2p^P3/2 — > 
ls^2p^Pi/2 transition for which the wavelength lies in the near infrared (IR) region, whereas the other two are in the IR 
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Table 2: Upper (lower) states /(i), transition wavelengths (A/i) in A, allowed (electric dipole) line-strengths (S/i), transition 
probabilities {Afi in sec~^) and oscillator strengths (f/i) in O VI. Numbers (x) in parentheses represent 10"^. RCCM labels 
the present calculation. The present calculation agrees well with the MCHF calculations done by Froese Fischer et al [8]. 
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JNIS 1 


GKASF 


/^TT ro 

L1V3 






NIST 


GRASP 


CIV3 


2 


1 


1037.63 


0.4502 


4.0824(8) 


4.09(8) 


4.2487(8) 


4.33(8) 


6.58961 


-2) 


6.60( 


-2) 


6.7338( 


-2) 


6.77( 


-2) 


3 


■J 
1 


1031.92 


0.9014 


4.1549(8) 


4.16(8) 


4.3212(8) 


4.39(8) 


A on/? cr / 

1.32651 


-1) 


1.33( 


-1) 


1.3552( 


-1) 


1.36( 


-1) 


A 

4 


I 


183.94 


A no ^ o 

0.0348 


5.6640(9) 


5.70(9) 


c o o /n\ 

5.6283(9) 


5.65(9) 


n o^nn/ 

2.87291 


n^ 

-2) 


2.89( 


-2) 


2.8726( 


-2) 


2.89( 


-2) 


4 


3 


184.12 


0.0700 


1.1361(10) 


1.14(10) 


1.1280(10) 


1.13(10) 


n O O T A / 

2.8870( 


-2) 


2.89( 


-2) 


2.8841( 


-2) 


2.90( 


-2) 


5 


1 


150.12 


0.0861 


2.5796(10) 


2.62(10) 


2.5791(10) 


O A / "1 A\ 

2.60(10) 


8.7158( 


-2) 


8.74( 


-2) 


8.7299( 


-2) 


8.79( 


-2) 


5 


4 


3835.30 


2.8158 


5.0564(7) 


5.05(7) 


5.1726(7) 


5.14(7) 


1.11501 


-1) 


1.11( 


-1) 


1.1224( 


-1) 


1.12( 


-1) 


n 

D 


1 


150.09 


0.1717 


2.5717(10) 


o oc\( A n\ 

2.62(10) 


2.5716(10) 


O /? A / 1 A\ 

2.60(10) 


1.73701 


-1) 


1.77( 


-1) 


1.7401( 


-1) 


1.75( 


-1) 


c 
u 


4 


QQ1 ceo 


t).U04/ 


o.ioUo^ / ) 




0.ZDD4(^ I ) 


O.ZO^ 1 ) 


Z.Z44DI 




2.24( 


-1) 


2.2590( 


-1) 


2.25( 


-1) 


7 


2 


172.93 


0.7451 


7.2975(10) 


7.33(10) 


7.3084(10) 


7.29(10) 


6.5437( 


-1) 


6.58( 


-1) 


6.5950( 


-1) 


6.58( 


-1) 


7 


3 


173.09 


0.1493 


1.4582(10) 


1.46(10) 


1.4604(10) 


1.46(10) 


6.5501( 


-2) 


6.57( 


-2) 


6.6009( 


-2) 


6.59( 


-2) 


7 


5 


11746.17 


3.8098 


1.1907(6) 


1.19(6) 


1.0478(6) 


1.01(6) 


4.9259( 


-2) 


4.90( 


-2) 


4.7036( 


-2) 


4.64( 


-2) 


7 


6 


11965.18 


0.7618 


2.2252(5) 


2.24(5) 


1.9794(5) 


1.91(5) 


4.8347( 


-3) 


4.81( 


-3) 


4.6144( 


-3) 


4.55( 


-3) 


8 


3 


172.99 


1.3439 


8.7667(10) 


8.78(10) 


8.7623(10) 


8.74(10) 


5.8994( 


-1) 


5.91( 


-1) 


5.9400( 


-1) 


5.93( 


-1) 


8 


6 


11894.88 


6.8581 


1.5331(6) 


1.37(6) 


1.2087(6) 


1.17(6) 


4.5388( 


-2) 


4.36( 


-2) 


4.1782( 


-2) 


4.12( 


-2) 


9 


2 


132.22 


0.0051 


2.2174(9) 


2.18(9) 


2.1448(9) 


2.18(9) 


5.8116( 


-3) 


5.70( 


-3) 


5.6488( 


-3) 


5.75( 


-3) 


9 


3 


132.31 


0.0102 


4.4464(9) 


4.34(9) 


4.2946(9) 


4.40(9) 


5.8349( 


-3) 


5.70( 


-3) 


5.6633( 


-3) 


5.80( 


-3) 


9 


5 


535.96 


0.2289 


1.5059(9) 


1.48(9) 


1.4819(9) 


1.51(9) 


6.4851( 


-2) 


6.39( 


-2) 


6.4073( 


-2) 


6.50( 


-2) 


9 


6 


536.41 


0.4725 


3.1011(9) 


2.96(9) 


2.9700(9) 


3.02(9) 


6.6885( 


-2) 


6.39( 


-2) 


6.4312( 


-2) 


6.53( 


-2) 



(ls^3p^P3/2 — > ls^3p^Fi/2 ) and far from the IR region (15^3^^1)5/2 — > ls^3(i ^£>3/2). These values are given in tables[3] 
andlH 

We have shown how the electric quadrupole {E2) and magnetic dipole (Ml) transition probabilities depend on the 
transition wavelengths in figure [3] and HI The abscissas in those two figures represent the different transitions involved. For 
example 3 — > 2 means transition from level 3 (ls^2p^P3/2) to level 2 (ls^2p ^Fi/2). The corresponding wavelengths are 
given in table [3] and [H On the other hand the ordinates in figure [3] and |4] correspond to the logarithms of A and f values 
respectively. From these two figures we have observed that for transitions between two fine structure levels, the Ml transition 
probabilities and hence the oscillator strengths are greater than those of the corresponding E2 transitions. Whereas in the 
other cases the trend is the opposite. This feature can be analyzed by an order of magnitude estimation following Landau 
and Lifshitz [17]. as given in the next subsection l|3.ip . 

Table [6] presents the comprehensive values of the lifetimes which are compared with the available theoretical data deter- 
mined by using the non-relativistic multi-configuration Hartree-Fock (MCHF) calculations by Charlotte F. Fischer's group 
[8] . We have calculated all possible decay rates of the excited states presented in figure [TJ Our calculated lifetimes are in 
good agreement with Froese Fischer's data except for the ls^4s^S'i/2 state which is dominated by the electric dipole allowed 
{E\) decay to the ls'^2p'^P^/2 state (given by 9 — > 3 decay in tabled]). In the present calculation we have considered all 
possible excitations of the core electrons to the virtual apace. This leads to a complete treatment of correlation effects in a 
more elegant and compact way. Higher order excitations are also taken into account through perturbative triple excitations 
which we have applied earlier and can be found in Refs. [16] llOj. Inclusion of higher order multipoles leads to the most 
accurate transition probability data sets of VI to our knowledge. 

Branching ratio of any decay of a particular state depends on the decay probability of that state to all the lower states 
and is defined as 

p _ ^up >low 

i up >low — .^-^ • (.J-^J 

/^^up ^low 

low 

We have analyzed the branching ratio of the 15^45^5*1/2 state to different lower states which decay dominantly via electric 

dipole {El) transition. The transition rates are given in table [2l Table [7] contains our calculated value of Tup >iow and 

the comparison with other available data. In the present calculation, the denominator in Eq. lfT2|) consists of the decay 
rates of the upper state (ls^4s ^S'1/2) to the lower states via all possible allowed and forbidden channels, whereas in the 
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Table 3: Upper (lower) states f{i), transition wavelengths (A/,) in A, electric quadrupole (E2) line-strengths (S/,), transition 
probabilities (A/^ in sec^^) and oscillator strengths (f/i) in O VI. Numbers (x) in parentheses represent 10^. 



/ 


i 


Xfi 


Sfi 


A/i 




3 


2 


187503.70 


0.4838 


5.8446(-10) 


6.1609(-15) 


5 


3 


175.68 


0.3540 


1.1845(6) 


2.7404(-6) 


6 


2 


175.47 


0.3519 


5.9233(5) 


5.4683(-6) 


6 


3 


175.63 


0.3615 


6.0557(5) 


2.8005(-6) 


6 


5 


641738.03 


18.6557 


4.7989(-ll) 


5.9256(-15) 


7 


1 


148.23 


1.0537 


4.1222(6) 


2.7157(-5) 


7 


4 


2891.26 


12.0391 


1.6683(1) 


4.1814(-8) 


8 


1 


148.15 


1.5807 


4.1336(6) 


4.0805(-5) 


8 


4 


2861.66 


18.0659 


1.7571(1) 


6.4714(-8) 


8 


7 


279530.06 


4.2353 


4.6320(-10) 


8.1388(-15) 


9 


7 


561.58 


2.2936 


2.3996(4) 


5.6725(-7) 


9 


8 


562.71 


3.5945 


3.5674(4) 


5.6449(-7) 



Table 4: Upper (lower) states f{i), transition wavelengths (A/j) in A, magnetic dipole (Ml) line-strengths (S/,), transition 
probabilities (A/j in sec~^) and oscillator strengths (f/j) in O VI. Numbers(a;) in parentheses represent 10^. 



/ 


i 


A/i 




A/i 




3 


2 


187503.70 


1.3333 


1.3639(-3) 


1.4378(-8) 


4 


1 


156.24 


5.808(-6) 


2.0538(1) 


7.5166(-11) 


5 


2 


175.52 


2.1609(-6) 


5.3898 


2.4894(-ll) 


5 


3 


175.68 


6.2410(-7) 


1.5523 


3.5915(-12) 


6 


2 


175.47 


2.9160(-7) 


3.6396(-l) 


3.3602(-12) 


6 


3 


175.63 


7.6038(-5) 


9.4641(1) 


4.3769(-10) 


6 


5 


641738.03 


1.3377 


3.4133(-5) 


4.2148(-9) 


7 


1 


148.23 


1.9600(-8) 


4.0581(-2) 


2.6736(-13) 


7 


4 


2891.26 


4.0000(-10) 


1.1160(-7) 


2.7974(-16) 


8 


7 


279530.06 


2.3997 


4.9393(-4) 


8.6792(-9) 


9 


1 


117.27 


1.4400(-6) 


1.2041(1) 


2.4828(-ll) 


9 


4 


470.24 


5.5225(-6) 


7.1627(-1) 


2.3746(-ll) 


9 


7 


561.58 


I.OOOO(-IO) 


7.6151(-6) 


1.8003(-16) 



Table 5: Upper (lower) states f{i), transition wavelengths (A/i) in A, magnetic quadrupole (M2) line-strengths (S/,), 
transition probabilities (A/j in sec~^) and oscillator strengths (f/i) in VI. Numbers (a;) in parentheses represent 10^. 



/ 


i 


A/i 


Sfi 


A/j 






3 


1 


1031.92 


6.8956(-3) 


2.1966(-5) 


7.0158( 


-15) 


4 


3 


184.12 


2.687(-2) 


9.4686(-l) 


2.4068( 


-12) 


6 


1 


150.01 


1.9721(1) 


9.6513(2) 


6.5210( 


-9) 


6 


4 


3812.52 


2.0035(1) 


9.2715(-5) 


4.0420( 


-13) 


7 


2 


172.93 


2.7019(-3) 


6.5114(-2) 


5.8407( 


-13) 


7 


3 


173.09 


8.3694(-4) 


2.0077(-2) 


9.0210( 


-14) 


7 


5 


11746.17 


3.2436(-4) 


5.4071(-12) 


2.2376( 


-19) 


7 


6 


11965.18 


1.7108(-4) 


2.6004(-12) 


5.5830( 


-20) 


8 


2 


172.83 


2.4806(-4) 


3.9977(-3) 


5.3723( 


-14) 


8 


3 


172.99 


5.4612(-o) 


8.7607(-4) 


5.8973( 


-15) 


8 


6 


11474.04 


3.8025(-6) 


4.7513(-14) 


1.4071( 


-21) 


9 


3 


132.31 


8.4971(-3) 


1.5621(0) 


2.0506( 


-12) 


9 


6 


536.41 


2.0101(1) 


3.3745(0) 


7.2804( 


-11) 
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Table 6: Lifetimes of VI states. MCHF labels Multi-configuration Hartree-Fock calculations [8|. Numbers (a;) in parentheses 
represent 10^. 



States 



IsVD 



3/2 



ISVD5/2 



Lifetimes (in sec 



RCCM 
2.44 95(-9) 
2.40 68(-9) 

5.87 37(-ll) 
3.86 88(-ll) 

3.88 05(-ll) 
1.14 20{-ll) 
1.14 06(-ll) 

11.04 55(-ll) 



MCHF 
2.44 60(-9) 
2.40 43(-9) 
5.85 46(-ll) 
3.78 54(-ll) 

3.81 14(-11) 
1.14 09(-ll) 
1.14 23(-ll) 
9.06 23(-ll) 



Table 7: Branching ratios of {up)ls'^4:s'^Si/2 — > {low)ls'^2p'^P3/2 , {up)ls'^As'^ S1/2 — > {\ow)ls'^2p'^Pi/2, 
(up)ls^4s ^5i/2 — > (low)ls^3p^P3/2 and (up)ls^4s ^5*1/2 — > (low)ls^3p ^Pi/2 transitions. The decays rates of the 
ls^4s ^5*1/2 states are dominated by electric dipole (El) transition. Each row corresponds to the value of the branching 
ratios of the upper level correspond to different El transitions obtained by different methods. The values given under NIST 
are obtained from the NIST database flSj. 



Methods 








^ls^4:S^Si/2 >ls^3p^Pl/2 


RCCM 


0.394 548 


0.196 737 


0.275 143 


0.133 610 


MCHF 


0.398 581 


0.199 117 


0.267 757 


0.134 543 


GRASP 


0.394 314 


0.196 928 


0.272 695 


0.136 063 


CIV3 


0.396 040 


0.196 220 


0.271 827 


0.135 914 


NIST 


0.395 985 


0.198 905 


0.270 073 


0.135 036 



other theoretical calculations only allowed transitions have been considered. To show the accuracy of our calculation we 
have presented the results upto sixth decimal places. Our calculated values of the branching ratios are off by 0.3% - 1.9% as 
compared to the observed value [18] . Among these transitions, the particular decay ((up)ls^4s ^S'1/2 — > (low)ls^2p ^P3/2) 
is the dominant one which is within an error of 0.3%. Similarly the branching ratios of decay of the other states can be easily 
obtained from this work and will be very useful for astrophysical and related studies. 

3.1 Order of magnitude estimation 

Following Landau and Lifshitz [17] the magnetic moment of an atom is equal in order of magnitude to the Bohr magneton 
II ^ — with all the variables with their usual meaning. This differs by a factor of a, the fine structure constant, from the 
order of magnitude of the electric dipole moment d ~ eoo ~ /me. Since vjc ^ a., we have approximated /i ~ dvjc. Hence 
it follows that the probability of Ml radiation from the atom is about o? times less than the electric dipole (i51) radiation 
at the same frequency. Magnetic radiation is therefore important in practice only for forbidden transitions governed by the 
selection rules for the electric case. The ratio of E2 to that of Ml radiation in order of magnitude for a given transition 
frequency w is 

E2_ {ea^fu?!^ a^m^o;^ /AP^ 

Ml ^ ^2 ^ ^\e ) ^^^^ 

where Ai? is the change in energy in the transition involved and the order of magnitude of the quadrupole moment and the 
energy of the atom [E] are ~ eag and ~ Y? jma^ respectively. For medium atomic transition frequencies (Ai? ~ E) the 
probabilities of E2 and Ml are of the same order of magnitude. If, however, Ai5 <C E, i.e. for transitions between two 
fine structure components of the same term the Ml transition probability is more probable than that of E2. The results 
presented in table [3] and [4] and in the figures [3] and [4] follow the pattern. 



7 



% error 




23456789 

States 

Figure 2: % error in calculation of relative transition energies of O VI. 




9^7 8^7 7^4 7^1 6 -> 5 6^3 6^2 5^3 3^2 



Figure 3: Dependencies of transition probability (A) on energies. The abscissa represent the energies for different transitions 
and the ordinate is the values of the A values in a log scale. Here 3 — > 2 means transition from level 3 (ls^2p^P3/2) to level 
2 {ls^2p^Pi/2). The levels 1 - 9 are identified as 1522525-^^2, ls^2pPi/2, ls^2p^P3/2, ls^3s^Si/2, ls^3p'^Pi/2, ls^3p^P3/2, 
1523^21)3^2, ls'^3d'^D^/2 and 1524525-^^2 respectively as denoted in table[l]in the column 'Index'. 
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4 Conclusion 

In this paper we have performed relativistic coupled-cluster calculations of all possible allowed and forbidden radiative 
transition probabilities of the first eight excited states of VI to the ground state. The lifetimes of those excited states have 
been determined considering all possible decay channels. This paper reports the first and most comprehensive calculation 
of the forbidden transition probabilities to our knowledge. Comparison with other available data are made for allowed 
transitions and radiative lifetimes. Accuracy of the computed transition data depends on both the accuracy of the state 
energies (transition energies) and the line strengths. Hence the determination of the state energies with an error of (and 
less than) 2% are used to benchmark the RCCM theory. We also made graphical analysis of the energy dependencies of the 
E2 and Ml transition probabilities and discussed the trends of different transitions with order of magnitude estimates. We 
believe that the present data sets will be useful for astrophysical studies of stellar spectra. 
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